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THE RADIAL SPANNING TREE OF A POISSON 
POINT PROCESS^ 

By Francois Baccelli and Charles Bordenave 
Ecole Normale Superieure and INRIA 

We analyze a class of spatial random spanning trees built on a 
realization of a homogeneous Poisson point process of the plane. This 
tree has a simple radial structure with the origin as its root. 

We first use stochastic geometry arguments to analyze local func- 
tionals of the random tree such as the distribution of the length of the 
edges or the mean degree of the vertices. Far away from the origin, 
these local properties are shown to be close to those of a variant of 
the directed spanning tree introduced by Bhatt and Roy. 

We then use the theory of continuous state space Markov chains 
to analyze some nonlocal properties of the tree, such as the shape 
and structure of its semi-infinite paths or the shape of the set of its 
vertices less than k generations away from the origin. 

This class of spanning trees has applications in many fields and, 
in particular, in communications. 

1. Introduction. There is a current interest in graphs generated over a 
random point set. This class of graphs includes, in particular, the minimal 
spanning tree, the nearest neighbor graph or the geometric graphs [20, 30]. 
On these graphs, the theory has mainly focused on large sample asymptotics 
when the random point set consists of independent uniform points in the 
unit square. 

In this paper we define a new model of rooted tree generated over a 
point set, the radial spanning tree (RST). This tree holds similarities with 
the minimal directed spanning tree [6, 21, 22] and the Poisson forest (as 
defined in [12, 13]). The mathematical analysis of the RST is conducted on 
a homogeneous Poisson point process (PPP) on the plane. The scope of this 
paper is two-fold: to prove some local geometric properties using stochastic 
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geometry and to derive global asymptotic properties. To this end, three 
embedded layers of analysis will be considered: local tree functionals around 
a vertex (e.g., the degree), properties of the path in the tree from a vertex to 
the root (e.g., the total length) or properties of the tree structure or topology 
(e.g., the characterization of the semi-infinite paths). These three layers of 
analysis rely on different probabilistic tools. 

This class of spanning trees has applications in many fields and especially 
in communication networks. For instance, in wireless sensor networks, the 
sensed information is gathered at a central point called cluster head and flows 
from the sensors along the edges of a spanning tree, the root of which is the 
cluster head [10]. Similarly, in multicast communication, the broadcasting 
of information flows from the source to the members of the multicast group 
along the edges of a spanning tree rooted at the source [5]. In both cases, 
the nodes on which the spanning tree is built are randomly located in some 
domain (sensors may be randomly spread in a region and may then move in 
a random way; similarly, the members of the multicast group are a random 
subset of the set of users in a large network). In both cases, it is crucial that 
the decisions to build the spanning tree from the random set of nodes be 
local; in addition, the properties of the path between a leaf of the random 
spanning tree and its root determine the network performance. Another 
application with the same characteristics can be found within the context of 
a class of information location algorithms used in peer to peer networks [7]. 

In the next section we give the basic definition and we summarize our 
main results. Section 3 gives various distributions on local properties on the 
tree. In Section 4 we focus on the directed spanning forest (DSF) which can 
be seen as the limit of the RST far away from the origin. In Sections 5.1 and 
5.3 we prove asymptotic shape theorems on the RST. Section 5.2 contains 
a proof of a law of large numbers on a semi-infinite path of the RST and 
Section 6 the proof of a law of large numbers for the spatial averages. Finally, 
in Section 7, we quickly discuss some extensions of the RST and give some 
open questions. 

2. The radial spanning tree. 

2.1. Definition. Let | • | denote a norm on and B{X,r) the open ball 
of radius r and center X. A set of points N of is said to be nonequidistant 
if there do not exist points X,Y,Z,T of N such that {X,Y} / {Z,T} and 
\X-Y\ = \Z-T\. 

If is a countable set of points in with no accumulation points, we 
write, for all bounded sets A, 

N{A):= ^ l{XeA). 
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Let be a countable set of points in W^, nonequidistant, with no accu- 
mulation points and such that O £ N. We define the RST of N,T = {N, E), 
with E the set of edges, as follows: each point, excluding the one in the 
origin, has an edge to its closest neighbor among the set of points which are 
closer to the origin. More formally, T is defined by 

If \Y\ < \X\ and X,Y £ N then {X,Y) G E 

^ N{B{0,\X\)nB{X,\X -Y\)) = 0. 

The nonequidistant property is needed to ensure that there is no tie: a 
vertex X which is not the origin has exactly one nearest neighbor which is 
closer to the origin. We trivially deduce that T is a tree. 

It is important to notice that the construction of the tree is local and that 
it does not minimize any global functional as do the minimal spanning tree 
or the trees analyzed by Howard and Newman in [17]. 

In this paper we will consider only point sets in the plane M^. All the 
results extend to higher dimension. If not otherwise specified, the norm | • | 
will be the Euclidean norm. We consider an orthonormal basis {0,ex,ey). 

Consider now some homogeneous Poisson point process N on the plane, 
with intensity A > and in its Palm version: almost surely (a.s.) O £ N. 
Since the Poisson point process is a.s. nonequidistant, we can a.s. generate 
the RST r of iV. 

An instance of such RST is given in Figure 1. 

Since T is scale- invariant, without loss of generality, we can set A = 1, 
for a general A, all results follow by multiplying distances by \/A- From the 
invariance of the PPP by rotation, we deduce also that the law of the RST 
is invariant by rotation. 

We now define another random graph which will be used in the analysis 
of the asymptotic behavior of the RST. Let (0,61,62) be an orthonormal 
basis of M^. On a nonequidistant point set with nonaccumulation, we 
define the directed spanning forest (DSF) with direction ei as follows: 
the ancestor of X G iV is the nearest point of N which has a strictly larger 
ei-coordinate. On a lattice, this random graph was first introduced by Bhatt 
and Roy in [6] and it also holds some similarities with the Poisson forest by 
Ferrari, Landim and Thorisson [12]. 

2.2. Notation. (A glossary of the main mathematical notation can be 
found at the end of the paper). The cardinality of set S will also be denoted 
by \S\. For all points X G N\0, the ancestor of X in the tree is the nearest 
neighbor which is closer to the origin. It will be denoted by A{X) and, 
by convention, we will take A{0) = O. The iterate of order A; of ^ will be 
denoted by A''{X). 

We will denote by Ro{X) the path from X to the origin, defined as the 
sequence of ancestors of X. In Section 5.1 Ro{X) will also be thought of as 
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a piecewise linear curve in R^, namely, as the union of the edges connecting 
the points of this sequence. 

Throughout the paper we will focus on functionals F defined on the ver- 
tices of T, such as the length of the edge {X,A{X)), its orientation, and so 
on. For any X G M^, F{X) is defined on the tree T constructed over iVU {X^. 
This definition is consistent with Slyvniak's theorem: if P is the probability 
measure of the Poisson point process A^, the Palm measure of A'^ with two 
points O and X is P * 6x * So- Hence, a.s. U {X} is nonequidistant and 
F[X) can be interpreted as the value of the functional F conditioned on the 
fact that X is a vertex of the tree. 

Several qualitative results of the present paper involve constants. For the 
sake of clarity, we will use Co to denote a positive constant to be thought 
of as small and Ci to denote a positive constant to be thought of as large. 
The exact value of Co and Ci may change from one line to the other and we 
could, for example, write Cq/Ci = Cq. The important point is that Co and 




Fig. 1. Radial spanning tree of 25000 points uniformly and independently distributed in 
the unit square. 
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Ci are universal constants that will never depend on the parameters of the 
problem. 

2.3. Summary of results. 

2.3.1. Topology of the radial spanning tree. We have seen that if is a 
Poisson point process under its Palm version, T is a.s. a tree. We will check 
also that this tree is a.s. locally finite (i.e., no vertex has an infinite degree). 

The next step to understand the intrinsic structure of T is to characterize 
its ends. An end is a semi-infinite self-avoiding path in T, starting from the 
origin (O = Yq,Yi, . . .). The set of ends of a tree is the set of distinct semi- 
infinite, self-avoiding paths (two semi-infinite paths are nondistinct if they 
share an infinite sub-path). 

A semi-infinite path (0 = lo> ^i) • • •) has an asymptotic direction if l^/|y„| 
has an a.s. limit in the unit sphere S^. 

The following theorem will be a consequence of Proposition 2.8 in [17]; it 
characterizes the ends of the RST: 

Theorem 2.1. The following set of properties holds almost surely: 

- every semi-infinite path has an asymptotic direction, 

- for every u £ , there exists at least one semi-infinite path with asymp- 
totic direction u, 

- the set of u^s of such that there is more than one semi-infinite path 
with asymptotic direction u is dense in . 

This theorem shows that the RST strongly differs from the minimal span- 
ning tree. In dimension two it has been proved that the minimal spanning 
tree has only one end (see [2]). 

Another property of the tree of interest to us is the set of points at tree- 
distance less than or equal to k from the origin T(/c) := {X E N ■.Af'{X) = 



Theorem 2.2. There exists a constant p> such that a.s., for all 
e > 0, for k large enough. 



O}. 





Moreover, a.s. and in L 



(2) 



In other words, the graph-diameter of the RST generated by a Poisson 
point process inside a ball grows linearly with the number of points. 
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2.3.2. Geometry of the radial spanning tree. In this paper we will focus 
on two types of geometrical results. 

The first type concerns functionals F{X,T) which depend only on vertices 
around X . For example, we will give explicit formulae for the distributions 
of L[X) := \X — A{X)\ (the length of the edge linking X to its ancestor), 
P{X) := \X\ — \A{X)\ (the progress to the origin) and the expectations of 
the degree D{X) of vertex X in the tree. 

More generally, we will prove a limit theorem for a large class of function- 
als, called stabilizing functionals. This class was first introduced by [18] and 
it was used by Penrose and Yukich [23, 24]; it is slightly modified here to 
suit to our framework. Roughly speaking, F{X,T) stabilizes T if the value 
at X depends only a small number of vertices around X . 

Since T depends only on the point set A^, with an abuse of notation, we 
can write equivalently F{X,N) for F{X,T) to make the dependence on N 
explicit. 

Definition 2.3. Let F{X,N) be a measurable M+-valued function de- 
fined on the vertices of a random graph G = {N, E) . 

F stabilizes G if for all X there exists an a.s. finite random variable 
R{X) > such that F{X, N n B{X, R{X)) U A^') = F{X, N) for any locally 
finite point set iV' C \ B{X, R{X)) and the distribution of R{X) does not 
depend on X. 

We will prove the following theorem where T-^^: denotes the DSF with 
direction —ex- 

Theorem 2.4. Let F be a stabilizing functional for T^e^- ^ tends 
to infinity, the distribution of F{xex,T) converges in total variation toward 
the distribution of F{0,T--e^) . 

Note that the subscript X in ex (the unit vector of the horizontal axis) 
has nothing to do with the real number x that we let go to oo. 

This theorem has to be related to the convergence of geometric graphs as 
it is defined for the objective method (refer to [1]). 

The second class of geometrical results is of a different nature: it concerns 
the path Ro{X) from X to the origin. The simplest result bears on H{X) := 
viii{k:A^{X) = O}, the generation of X in the RST, that is, the cardinal 
of Ro{X). Let Xq = X,. . . ,Xff^x) = O be the sequence of points of N in 

Ro{X). Along the line of [17], it is interesting to look at X)fc=o l-^fc+i — 
with a > 0. More generally, we will prove the following: 
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Theorem 2.5. Let p be the constant in Theorem 2.2. There exists a 
probability measure n onM. such that if g{X) is a measurable function from 
to M, \g{X)\ < max(Ci, for some a>0, then a.s. 

HiX) 1 1 ^^-^^ 

ix^^ioo ^ = p H{X) g ^^[^^-^ - ^'^l^-^) = 

where \U]v is the vector U rotated by an angle —0 and V = rcosO ■ ex + 
rsin^-Cy is the vector u expressed in the local coordinates ofV). 

If g is continuous, we also have 

|X|-+oo H[A) 

We prove in Section 4 that the probability measure vr can be interpreted as 
the stationary measure on the infinite edge process in the DSF. Theorem 2.5 
is a law of large numbers and vr can be understood as the limit probability 
measure of an edge conditioned on being on a semi-infinite path. 

Theorems 2.4 and 2.5 are of a different nature. In particular, the mean 
value of L{X) differs from the average of the lengths of the edges along the 
path Ro{X). This can be explained by the fact that being on a long path 
is a bias. We will discuss this in Section 6. 

3. Local properties of the radial spanning tree. 

3.1. Distribution of the length of edges. Let X G and L{X) := \X — 
A{X)\. Let < r < \X\; we get 

¥{L{X) >r) = l(r < \X\)¥{N{B{X,r) n B{0, \X\)) = 0) 

(3) 

= l(r< |X|)e-^'^(l^l'^), 

where M(x,r) is the volume of the lens of the right part of Figure 2. Using 
the formula for the surface depicted by the left figure of Figure 2, we get 
that 

(4) M(x, r) = xi<P- ^^2f^__t_^j^ 



-2 2 
with 

r 

6 = 2 arcsin — . 
^ 2x 

With an abuse of notation, for x > 0, we define: L{x) := L{{x,0)). L{x) 
has a density on (0, x) equal to 

dr 
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and a mass at x equal to 



(5) 



-M{x,x) 



-a:2(27r/3-sin(27r/3)) 



2{27r/3-V3/2) 



Notice that the distribution function of L{x) is not stochasticahy mono- 
tone in X. Its mean E(L(x)), which is plotted in Figure 3, is not monotone 
in X either. 

Remark 3.1. In Appendix, Section A.l, we also consider the distribu- 
tion of L{Xn) when we number the points of N by their a.s. increasing 
radius. 

3.2. Distribution of progress. Given L{X) = r < \X\, consider the angle 

eiX) := OXMX). Since ip = 7r/2-0/2 (see the right part of Figure 2), e{X) 
is uniformly distributed on the interval (tt -|- arg(X) — ^,7r + arg(X) + 
with cost/; = sin(i?!)/2) = r/{2\X\), that is, t(j = arccos 2\x]- Given L{X) = \X\, 
the angle 0{X) is vr -|- arg(X). 

The joint distribution density of {L{X),6{X)) is equal to 



(6) 



l{r e (0, |X|))^M(|X|,r)e-^-^(l^l'") dr 
dr 

X 1(6' G (vr + arg(X) - V', vr + arg(X) + -0)) 



d£_ 

24) 



~M(\Xl\X\) 



+ ^\X\{f)5lT+^Tg{X){ff)(i' 

The progress is defined as 

P{X):=\X\-\A{X)\. 
The mean progress P((x,0)) is plotted in function of x in Figure 3. 




Fig. 2. Left: the surface of the dashed lens is equal to ■'-^ |26 — sin26| . Right: the dashed 
lens. 
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Fig. 3. Left: Mean of L{x) in function of x. Right: Mean of P{x) in function ofx. 



3.3. Mean degree of a vertex. 

3.3.1. Degree at the origin. The degree of the origin is 

D{0)= t{N{B{T,\T\)r^B{0,\T\)) = ^). 

T&N\{0} 

Hence, using Campbell's formula, we get 

(7) ED(0) = 2tt / e-'- (2V3-sm{2V3))^ ^ ^ 2.56. 

Jo 27r/3 — v 3/2 

The following property is also of interest. In the same vein as the bound- 
edness of the degree of a vertex in the minimal spanning tree, we have the 
following: 

Lemma 3.2. The degree of vertex O is hounded from above by 5 a.s. 

Proof. Order the points directly attached to the origin by increasing 
the polar angle. Let X and Y denote two neighboring points in this sequence. 
Assume \OX\ < \0Y\. Denote by (p the angle between these two vectors. We 
have 

ixVp = \OX\^ + |Oy|2 - 210X110^1 coscj). 

Since Y is attached to the origin, necessarily |Xyp > |Oyp, which implies 
that 

2|ox||oy|cos(?;>< |oxp. 

Using now the assumption that \0X\ < \OY\, we get coscp < 1/2. Hence, 
|(/>|>7r/3. □ 
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3.3.2. Degree outside the origin. The degree of a vertex X ^ O is given 

by 

D{X) = 1+ H\T\>\X\)1{N{B{T,\X -T\)nB{O,\T\)) = 0) 

T£N\{X} 

xl{0(^B{T,\X -T\)). 

Indeed, a point T of modulus larger than \X\ shares an edge with X if and 
only if there is no point of smaller modulus closer from T than X. 

Let X ^ O, \X\ = X and denoting ¥,D{X) by ¥,D{x), using Campbell's 
formula while taking the expectation of equation (8), 

ED{x) = 1 + E ^ 1{N{B{T, \X - T\) n B{0, \T\)) = 0) 



X i(x< |r|)i(|r| > \x-T\) 

f fa.vccos{x/{2p)) 

1+ / / e^^^^'^P^^^pdpdO, 

J p>x J ~ arccos(x/f2p)) 



farccos(x/ (2p)) 
I p>x J ~ arccos(x/(2p)) 

where Q{x,p,9) is the dashed surface in Figure 4 for X = (x,0) and T ■ 



{p,0). The condition that \T\ > \X — T\ [or, equivalently, that 9 belongs 
to the interval (— arccos(^), arccos(^))] translates the fact that the origin 
should not be contained in this lens. Hence, 



/•arccos(a;/(2p)) „ , , 
ED{X) = 1+ I / g-pV2|2a-sin2a| 



p>x J — arccos(a;/(2p)) 

^ g-(p2+a;2-2pa;cos0)/2|2/3-sin2/3|^^^^^^ 

where a and /3 are the angles depicted in Figure 4. 

If M = ^, we have cos a = (1 — m~^)/2 + cos 6 and (3 = {tt — a)/2. Fi- 
nally, 

/■ /■arccos(l/(2u)) „ „ 

ED(X) = l + 2xU / ^-{u^x^)/2i2a-sin2a) 



(9) 

^ ^-xy2{l+u^-2ucose){n-a-sina)^^^ _ 

The mean degree is plotted in Figure 5. 

The following lemma is remarkable in view of Lemma 3.2. 

Lemma 3.3. The degree of a vertex of the DSF is not hounded and in 
the RST a.s. 

sup D{X) = +00. 
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Fig. 4. Left: Q{t,p,9). Right: The a and [3 angles. 




Z 4 I ■ iq 

Fig. 5. ED(y) as function ofy>0. 



Proof. Let A-e^{X) be the ancestor of X in the DSF with direction 
-e^,;, T_e,. The DSF built on the point set {X„ = (2-",3''),n G N} U {0} 
gives, for ah n, A-e^{Xn) = O, in particular, the degree of the origin is 
infinite. 

We now prove the second statement of the lemma. 
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Let M G N*; for n > 0, we define U„ = [2"" - e, 2"" + e] x [S"" - e, 3" + 
e], U_i = [-e,e] x [-e,e], Am = B{OA^') \ (U_i<n<A/Un) and Em{X) = 
{N{X + Am) = 0,N{X+Vn) = 1, -1 < n < M}. We have P(-Ea/(X)) = 5>0 
and if \X-Y\> 2.4*'^, £'m(^) and £;m(>") are independent (for e is small 
enough) . 

For £ small enough, if Em{X) occurs, the point in U_i(X) has degree at 
least M in T-f^^. Similarly for the RST, if \X\ is large enough and if Em{X) 
occurs, the point in \]-i{X) has degree at least M in T. 

Using the independence of the events EM{'^kA^ ex), A; G N, we deduce that 
these events appear infinitely often and this concludes the proof. □ 



3.4. Limit distribution of an edge length. Using equations (3) and (4), a 
direct computation gives 



lim P(L(X) > r) = exp , 

lim mL{X)]= I e-^-'^/^dr- ^. 



By similar arguments, the asymptotic progress has for Laplace transform 

2^ POO /■7r/2 / ^ 

|X|^+oo ' ' TT Jr=oJe=~w/2 

in: 



lim E(e-^^W) = -/ / e-^'^^'^^exp )TTrdrd9 

H+oo TT Jr=oJe=~TT/2 V 2 / 

/•CO /•7r/2 / 7-j,2 \ 

/ / e-^^'^^^exp )rdrde. 

Jr=0 Je=-Tr/2 V 2 / 



/r-=0 J0=-7r/2 

In particular, the mean asymptotic progress is 

V2 



(12) lim EP{X) = - / exp cos dr = — . 

\X\^+(xi vr Jr=0 Je=-7r/2 V 2 / TT 



L and P are distributed as the length of the edge {0,A{0)) and its x- 
coordinate in the DSF. This result is contained in Theorem 2.4 which we 
will prove in Section 3.6. With the same limit reasoning, for the degree, we 
get 

lim ED{X) = 1+ f expf-^^i^) (iy 

\X\^+oo J Half-plane V 2 / 

rTT/2 r+oo / ^^2^ 

= 1 + / / exp \rdrd6 

J-./2J0 V 2 J 

= 2. 
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3.5. Expectation of the number of crossing edges. Let x > and C{x) be 
the number of edges connecting a vertex belonging to the closure of i?(0, x) 
and a vertex outside. Almost surely, we can sort the points of N by increasing 
norm O = \Tq\ < |ri| < • • •. We have 

C{\Tn+l\) = C(\Tn\) + D{Tn+i) - 2. 

It follows that 

C{x) = D{0)+ J2 iD{Tk)-2). 

\n\<x 

Taking expectation, we deduce from Campbell's formula that 

EC{x)=ED{0)+2tt r {'KD{t)-2)tdt. 
Jo 

In Section 3.4, we proved that limx->+oo ED{x) = 2. Hence, the evaluation 

of the asymptotic behavior of EC(a;) r6C|iiirGS that of c — — liin^; >-i~oo 

We would then deduce limx^^oc^C (x) / x = c. The exact computation of c 
is beyond the scope of this paper. 

Let 7 be an arc on the circle of radius r and center O and let C(r, 7) be 
the counting measure of the number of edges crossing 7. From the invariance 
by rotations of the PPP, we deduce that EC(r, 7) = ^(7)^^^, where ^(7) is 
the length of the arc. In other words, the point process of edge crossings on 
the circle of radius r is stationary and with intensity /z(r) = ^^^J^^ . 

3.6. Proof of Theorem 2.4. Let X G \ {0} and A-eA^) be the an- 
cestor of X in Tlgj with direction — ei (a given vector). We define 9 as the 
angle between ei and X. The next lemma compares A{X) with A-ei{X) if 
we build T and Tlei on the same PPP. 

Lemma 3.4. Let X £R'^\ {O}; there exists C\ such that 

p(^(x) / A-eA^)) < 1 A Ci{i/\x\ + e). 

Proof. Without loss of generality, we suppose X = xcx, x > and 
9>0. The sets L{X,ei) and K{X,ei) are depicted in Figure 6. We note 
that iiAdiX)^K{X,ei) and AiX)^L{X,ei)U{0}, then A^ei{X)=A{X) 
and, hence, 

P(^_ei {X) / A{X)) < P(^_ei {X) G K{X, ei)) 

+ F{A{X) G L{X, ei)) + ¥{AiX) = O). 

The last term is easily computed: F{A{X) = 0) = e~*'^(^'^). 
To upper bound the second term, we notice that L{X,ei) is contained in 
a cone of angle 9 (see Figure 6); hence, W{A{X) G ei)) < 9/tt. 
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The first term is upper bounded similarly. Let K-{X,ei) be the lower 
part of K{X,ei); we have 

P(^_ei(X) G K{X,ei)) < 2P(^_ei(X) G K-iX,ei)) 

< J 2-Kre~'^'''^ ^arcsin^^^ + 0^ dr 

<Ci{i/x + e). □ 

We now prove Theorem 2.4. We recall that F{X,T) is the value of the 
functional F at vertex X when the RST T is built on NU{X}. Without 
loss of generality, we can suppose that X = xcx, x >0. We build T and Tle^ 
on the same PPP. For all r > 0, we define the event Jr{X) = {T n B{X, r) = 
T^ex n B{X,r)}. Let F be a stabilizing functional for T-^x- Using the ter- 
minology of Definition 2.3, 

P(F(X,T)/F(X,T_eJ) 

<njR[x){xT) 

<F{R{X) > r) + F{Jr{Xf) 
<nR>r)+n U A{T)^A-ex{T)\ 

\T(^Nr\B{x,r) / 

< F{R > r) + F{N{B{x,r) >n)+ nCi{r/{x -r) + l/{x- r))), 

where we have used Lemma 3.4. For r/ > 0, we fix r such that F[R > r) <rj. 
Note also that F{N{B{X,r)) > n) < exp(-nln^) (Lemma 11.1.1 of [28]). 




Fig. 6. The sets L{X,ei) and K{X,ei). 
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Hence, taking n = [x'^l for some < a < 1, we deduce that 
limsupP(F(X,T) / F(X,T_eJ) < v, 

and it follows lim|x|^ooIP'('^(^,'^) / F{X,r_eJ) = 0. 

To complete the proof, notice that 7^ is stationary along e-x ■ F(xe^, Tle^;) 
and F{0,T-e^) have the same distribution. 



Remark 3.5. It is easy to check that the vector (L{X),9{X)) [and hence 
the progress P{X)], or the degree D{X) are stabilizing functionals for T-g^. 
So are the first k segments of the path from X to the origin in the DSF, for 
all finite k, or the subtree of the DSF rooted in X and of depth k. 



4. The directed path. The directed spanning forest is the key tool for 
the asymptotic analysis of the RST. This section is dedicated to the analysis 
of the paths originating from some vertex. 



4.1. The underlying Markov chain. We denote by To,Ti,...,T„ the se- 
quence of the successive ancestors of point Tq in the DSF with respect to 
the direction —ex- For i>l, the edge between Tj and Tj_i is Ui := Tj_i — Tj. 
Let Uq := —Tq. From the definition of the tree, for i > 1 the x-coordinate of 
Ui is positive (almost-surely) and we have: 

n n 

Tn = -^Ui and Tn-Tk = - ^ Ui for < n. 

1=0 l=k+l 

The conditional distribution of Un+i given (C/q, . . . , Un) can be analytically 
determined. Indeed, let 7i be the half-plane (x > 0), Dq := 7i and for n > 1, 
let 

{n-l ~j c 

\jBiTr,-Tk,\Uk+i\)\ 
k=0 J 

(see Figure 9). Given {Uq, . . . ,Un), Un+i admits the following density in 
polar coordinates: 

(14) #|i?(O,r)nD.|e-im0nZ).| Hrmn^-drde 



dr' ' ' ' ui{C{0,r)nDny 

where z/i denotes the 1-dimensional Lebesgue measure and C{T,r) the circle 
of radius r and center T. 

Let Dn be the set given by equation (13) for some fixed points Tq, . . . , T^. 
The following two lemmas will have important applications. 
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Lemma 4.1. Let {rn,9n) denote the coordinates ofUn- Then for all < 
a < 7r/2, 

P(|0„| <a|D„_i,r„)>-. 

TT 

For —tt/2 < a < P < tt/2, if the cone {{r, 6); 9 £ [a, (3]} is included in Dn-i, 
then 

P(0„G[Q,/?]|D„_i,r„)>^^. 

vr 

Proof. Fix r > 0, in view of the geometry of Dn-^ 

[O e (-|, I) : (r, 9) G = (-0_(r), 0+(r)), 

with 0<6'± < f . 

From equation (14), for a fixed r, the p.d.f. of 9i conditioned on rn = r 
is uniform on (— 0_(r), 0+(r)). li 9- <a and 0+ < a, then P(|0n| < a^n = 
r,Dn-i) = 1; else, supposing, for example, 0+ > a, then 

P(|a„| < a|r„ = r, Z)„_i) > P(0 < 9n < a\rn = r, D^-i) 

a a 
> > -. 

- 9+ + 9_- TT 

The proof of the second assertion is similar. □ 

Define the cones 

^( X = {r,9)eR'^:9e[0,a)}, for a > 0, 
\X = {r,9)eM?:9e la,0], for a < 0. 

Lemma 4.2. For all n, 

C Dn or C_^/6 C Dn- 

In particular, if {rn,9n) denote the coordinates ofUn, then 
P(r„ > u\Tn.uTn-2, . . . ,To) < e--"'/i2_ 

Proof. The proof relies on a simple geometrical argument. Suppose first 
n = 2 and consider a circle C2 of radius r2 and center O (set to be Ti) and 
another circle Ci of radius ri and center Tq, O £ Ci. In polar coordinates, 
To is at {ri,9i) and T2 at {92, r 2). For 1^2 — 6*1! < f j the equation of Ci 
in polar coordinates is r = 2ricos(0 — 9i) (see Figure 7). The point T2 is 
somewhere on C2; suppose, for example, that it is in the orthant 92 G [— f , 0]. 
If 9i ^ [— vr,— ^] or r\ > 2r2Cos^i, then D2 contains c^/2- 
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Suppose instead 6i G [— vr, — ^] and ri < 2r2 cos^i (see Figure 7). We have 
to prove that the largest cone with origin T2 contained in D2 contains c^r/e 
or c_^/6. 

The worst case is when T2 is at M , defined as the intersection of Ci and C2 
in the orthant 6e [— 1,0]. We have M = (r2,(^), with = + arccos(^). 
In this case, an easy calculation shows that the largest cone contained in 
D2 with origin M = T2 is {X = M + (r, 9) : 9 £ {-^ - 9i + 2^, f + </>)}. 
The worst case is reached when = ^, since 9i > — tt. We deduce that 
max(f + ^i-2(/),f +(/.)> f. 

This concludes the proof for n = 2. For n > 3, the largest cone contained 
in Dn with origin T„ is tangent to (at most) two circles, and the same 
conclusion holds. □ 

A sample path together with the associated exclusion discs are given in 
Figure 8. 

From equation (14), the process {Un},n G N, is not Markov. We may 
circumvent this difficulty by defining tq := 1 and 

{r m-2 ^ c 

m>Tn:nnl U B{T^-Tk,\Uk+i\)\ =0 
U=r„-1 J 

In words, ri is the first time m such that no disc previous to m — 1 intersects 
the half plane at T^ . 

Each Tn is a stopping time with respect to the process history generated 
by {Um}- We call these times Markovian times (we will soon see why). 




Fig. 7. D2 and the largest cone with origin T2 contained in D2- 
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S3 m M MS U Mil 

Fig. 8. A sample path and the associated discs. 



Let Li := denote the length of the edge from Tj to Tj+i and Pi := 

{Ti — Ti+i,ex) the directed progress along this edge (where (•, •) denotes the 
usual scalar product). 

The projection of Ti.\ =Ti. H {B{Ti — Tq, on the x-axis is an 

interval of length S^i = Lq — Pq. More generally, denote by ^„ the projection 
of H\ Dn = H D {Ufc=d ^C^n — Tk, \Uk+i\)} on the x-axis. This sequence 
satisfies the recurrence relation 

(16) in+l = max(^„ - Pn, Ln - Pn), 

so that, for all n > 1, 



n> 1, 



max Lj_i 
l<i<n \ 



n-l N 
k=i—l / 



and the Markovian time ri is then simply rewritten as 

Ti = inf{m >l:Pm> Cm}- 
For instance, on the realization of Figure 9, ti = 4. 

Lemma 4.3. ForaUm>l, 
P(r^ - Tm-i > n\Dr,) < Cie-^f"^ and P(r™ > (3m\Dr,) < Cie 



Com 
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Fig. 9. The process To,Ti, . . . ,T„, the associated Ui,. . . ,U„ and in dashed H \ D2- 

for some positive constants P,Cq,Ci (depending on Dt-q), so that, in partic- 
ular, Tm is a.s. finite. 

Proof. This result follows from the statement of Lemma A. 2 (in Ap- 
pendix). We explain the connection between the setting of Lemma A. 2 and 
the current setting in the particular case m = l. 

We first show that, for all n, the random variable ^„ defined in (16) is 
bounded from above by where is the maximal residual service time 
just before the nth arrival in a GI / GI /oo queue with i.i.d. service times 
{Ln} and i.i.d. inter-arrival times {Pn} where: 

• Li is the distance from Tj to the closest point in the cone (either c^^/q 
and c_t^/q) that is fully included in Di if there is only one such cone. If 
both c^/g and c_^/g are included in Di, then one selects one of them at 
random; 

• Pi is defined as follows: one samples an independent Poisson point process 
A^' of intensity 1 and one picks the point of 

{{N n Di) u {N' n Df n n)) n BiTi, Li) 

which has the smallest progress from Tj. 

Since for all i, Li < Li and Pi > Pi a.s., an immediate induction gives < 
a.s., where ^i is defined by the same recursion but with Li and Pi in place of 
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Li and Pi. This recursion is that of the maximal residual service time just 
before the arrivals in a GI / GI /oo queue. 

Using a Loynes type argument (see, e.g., [4]), the sequence {.^j} is eas- 
ily seen to be stochastically monotone in i and to converge weakly to a 
nondegenerate limit when i tends to infinity. 

Pick a such that < a) = 77 > 0. Let {rn,0n) be the polar coordinates 
of Un ■ In view of the remark preceding the lemma, 

P(ti = m\Do) > P(em < a)P(r,„+i cos(0^+i) > a|D„) 
> P(|< a)P(r„+icos(0„,+i) > a\D^). 

Let c(r, a) denote the set {{p,6) : p < r;\9\ < a}. From Lemma 4.1, for all 
sets Dm, for all a £ (0, |), 

P(rm+i cos 6*^+1 > a\Dm) > F{rm+icos0m+i > a; \9n\ < a\Dm) 

> ^g"i^2(c(a/cosQ,a)nD,n) 

~ vr 
vr 

Let 

fi = inf{m >lPm> im}- 

Since ri < fi, in order to prove the finiteness of ti (or the exponential bound 
on the tail of its law), it is enough to prove this property on fi and this is 
Lemma A. 2. □ 

For n > 1 , we define the path between two successive Markovian times as 
the point process 

(17) = {0,Tr„_, - Tr„_,+i,. . . ,r,„_i - TrJ 

and we define <I>o = {0,To — Ti}. From Lemma 4.3, for all n, is a finite 
point process on ^ = x M. We endow the set of finite point processes 
on X with the usual weak topology, {J\fx-,B{J\fx))-, making a complete metric 
space (see [8] for details). 

Theorem 4.4. {$„},n G N, is a positive and uniformly ergodic Harris 
chain. 

The proof is organized in several steps. All definitions on Markov chains 
are taken from [19]. 
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Markov property. Almost surely, we may write = {0, , . . . , X^}, 
where the x-coordinates of the X" are increasing. By the definition of r^, 
for all n, we have 

Dr,,=nnBi-Ur^,\Urjr. 

It follows from equation (14) that the law of f/r„+i = depends on ^q, . . . , <I>, 
only through Ur,^ = T^„_i - T^^ . 

More generally, since H n B{Tn - T^, \Uk+i\) = imphes H n B{Tm - 
Tfc, |f7fc_|_i|) = for m>n, it follows that, for all m > Tn, 

{m-l ~| c 

y BiTm-n,\Uk+i\)\ , 
k = Tn-l ) 

which depends only on Ur„, ■ ■ ■ , Um- In particular, for all n, 

P($„+l G .|^0, . . . , ^n) = IP($n+l G Hf^rJ = F($„+i G .|^™), 

since U-r„ is a functional of only, which completes the proof of the Markov 
property. 

In what follows, it will be useful to consider more general initial conditions 
than the two-point $o defined above. Any finite point process <I> on X with 
$ = {0, Xi, . . . , and satisfying the constraint 

DN = 'Hr\B{-UN,\UN\Y 

is an acceptable initial condition when taking tq = N. A special case is that 
where $ = {0}, where = and Dq = 7i. 

In what follows, for a point pattern $ as above, we will also use the nota- 
tion Aj($) for the ith point (the point, are ordered by their nondecreasing 
x-coordinates), A^($) for the number of points and -D(<I>) for the set Dat. 

As commonly done in Markov chain theory, we will denote by P$ the 
probability P(»|<l>o = ^) and by the probability P(<I>„ G •|<I'o = $), with 
<I> = {0, Ai, . . . , X]\f} a finite point process satisfying the above property. In 
particular, Pq (or Pq) will denote the probability measure of the process 
conditioned on $o = {0}. 

Note that P$ depends on $ only through D($). 

Irreducibility and aperiodicity. For A £ B{J\fx), define 

(18) cr^ :=min{n> l:^>n G A} 
and 

L{^,A) := P{aA < oo|$o = ^)- 

The irreducibility of a Markov chain on J\fx relies on the existence of a 
measure on I3{Mx) such that, for all A in B{J\fx), 

(19) v{A) > implies L($, A) > for all <^>eAfx- 
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For all sub cr-algebras Q of B{J\fx), let P^lg denote the restriction of 

to g. 

Let denote the sub cr-algebra of B{J\fx) generated by the sets {N{(^) = 
l}n{Xi($) G e 0+, where 0+ denotes the positive orthant M+ xR+. 

We choose v to be -Pq^Ijf. 

As already pointed out, depends on <I> only through D{^). In view 
of equation (14), if D{(^) n 0+ = 0+, then P^|jc- and are equivalent 

measures. So if D{(^) 0+ = 0+ and Pq{A) > 0, for some ^ e JF, we then 
have P^iA) > 0, so that A) > 0. 

Consider now the case where D{^) f] 7^ O'^ . The first point of $1 
is Xi{^i) = Uro+i- Let (Pi, 61) be the coordinates of Xi{^i). Let (0,^) 
denote the projection on the x axis of the set 7i \ D{^). From equation (14), 
P(Pi cos 9i > ^, 61 > 0) > 0. If Pi cos 61 > ^, To + 1 is a Markovian time, so 
that n = To + 1, and D{^i) = nn B{-Uri,\Uri\y . If, in addition, 61 > 0, 
P'(<I>i) nO+ = 0+, from what precedes, for A G and such that Pq{A) > 0, 
then P|(A) > 0, which implies L{^,A) > 0. 

The irreducibility is thus proved. The proof of the aperiodicity is along 
the same lines. 

Small set. A set S G B{Mx) is small if there exists an integer n> and 
a nontrivial measure v on B{J\fx) such that, for all A G B{Mx), 

(20) mf P|(^)>K^). 

For $ = {0,Xi, . . . in Afx, we define (P(<I>), 6($)) as the coor- 

dinates of Xr - Xr-i- For r > and < a < 7r/2, let S{r,a) G B{Mx) be 
defined by 

(21) 5(r,a) ={^(^Mx: < a, P(^) < r}. 

Let Q denote the sub fi-algebra of B{J\fx) generated by the sets {r(<I>) = 
1} n {Xi{<^) e B},B eCa, where Ca denotes the cone 

Ca = {ir,6)eR'^:ee [-tt/2 + a,7r/2 - a]}. 

We now prove that 

(22) S = S{r,a) 

is a small set for v = g{r,a)PQ\g and n = 1, where g{r,a) is a positive 
constant to be determined below. 
From equation (14), for all <I> G 5, 

P$(Xi($i) G C«,t($i) = 1) = P$(|ei| < tt/2 - a. Pi cosGi > r) 

> ^'o(|0i| < vr/2 - a, Pi cos Gi > r) 
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Note also that from equation (14), for all B C £, 

P^{Xi£B\\ei\ <Tr/2-a,Ri cosGi >r) 

= Po{Xi{<Pi) e B\\ei\ <7r/2-a,i?icosGi>r). 
These two remarks imply that, for all A(^Q, 

P^i<Pi£A)>g{r,a)Poi<^i£A), 
which concludes the proof. 

Positivity and uniform ergodicity. Let as = inf{n > 1 : G 5} be the 
first return time to S with S defined as above. 

Lemma 4.5. For all r > 0, S = 5(r, 7r/6), we have sup^ E^ag < oo. 
More precisely, sup$P(cr5 > n) < (1 — 6{r))"' , with 6{r) = (1 — e~'^^ /^^)/6. 

Proof. From Lemmas 4.1 and 4.2, for all n, 



from which one easily deduces that sup,j, E^as < oo. □ 

In view of Theorems 10.4.10 and 16.0.2 of [19], is positive and uni- 

formly ergodic. 

4.2. Limit theorems. We follow the approach of Athreya and Ney (see 
[3]). Since the Markov chain $ is Harris recurrent and strongly aperiodic, we 
can build an increasing sequence of finite stopping times Nj. (on an enlarged 
probability space) such that A'^o = and 



where /i is a probability measure on Mx- Therefore, N^^ — 1 is a regenerative 
time: the sequences ($„_|_7Vfe), ?t. G N, and ($n)>0 < < A^^ — 1, are indepen- 
dent. 

Lemma 4.5 implies that — N^^i is stochastically dominated by a geo- 
metric law. We define 9^ = tat^ with defined in (17); since ^q^, < r [with 

defined in (16)], from Lemma 4.3, P(0fc+i — 6k> t\!FQ^) < Ci exp(— Cq/c), 
for some positive constants. 

Finally we have the following: 




Hence, for all ^, 



(23) 




(24) 



P$„ ($„ G A, Nk =n)= /i(A)P$o {Nk = n), 
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Theorem 4.6. There exists an increasing sequence of finite stopping 
times {6k}, k G N, (on an enlarged probability space) such that 9q = 0, for 
all k > 0, F{6k+i — Ok > n\Tgf,) < Ci exp(— Con), for some positive constants 
and the vectors {Ue^^i, . . . , Uei^^-^),k G N* are i.i.d. 

For ^ = {0,Xi,...,Xr} in Afx, we define /($) = ^^=1 niax(l, |X„|"), 
a> 1. Using the upper bounds of Lemma 4.3 and Lemma 4.2, we get 

(25) supExf E /($fc)) <oo, 



k=0 



where as is defined by equation (18) and S is the smah set defined above. 

Let n be the invariant distribution of the Markov chain. We define the 
invariant distribution of the edge process Tq — Ti, . . . , T„_i — r„, . . . as 

(26) 7T{A)=En{Tr^En(j2l{{Xi+i-Xi)eA)^, AeR^. 

We may now deduce a first limit theorem. 

Theorem 4.7. For all measurable functions g, let Sn = Sfc=o dC^k+i — 
Tk). If g{X) < max(C, \X\'') forC>0, a> 0, then a.s., 

]ira^=En{rr'EJj29iXi+i-Xi)] =7r{g). 



Let G($) = E[=o (9(^m - Xi) - <9)); if 



7^ = En{G{^of) + 2 En(G($o)G($fc)) > 0, 
fc=i 



then a central limit theorem also holds: 



^ (5„-^(5))^Ar(0,l). 



7Vn 

lf"f = 0, then -^{Sn — T^{g)) tends a.s. to 0. 

This theorem characterizes the asymptotic behavior of a path in the di- 
rected spanning forest. 



Proof. Theorem 4.7 is a direct application of Theorem 17.0.1 of [19]. 
Theorems 14.2.3, 14.2.4 and 14.3.7 of [19] along with equation (25) ensure 
that all requirements are fulfilled. □ 
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Corollary 4.8. There exists positive constants p, py and for a > 0, Ic 
such that 

2 n— 1 



(27) lim - V (Tfc - Tfc+i , Ca;) = = p, 

n— ♦CO 77, ^ — ' 

(28) \mi -y^\{Tk-Tk+i,ey)\=py, 



n — s-OO 77 



I n—1 

(29) lim -^irfc+i-T^r = /„. 

By simulation of 20000 transitions of the chain, one obtains that p ~ 0.504, 
Py ~ 0.46 and li ^ 0.75. The value of p is significantly larger than the mean 
asymptotic progress as evaluated in (12). The latter can be seen as the 
expectation of the progress from point Tq = under the Palm probability of 
the Poisson point process N , whereas p is the expectation of the progress 
under another measure: for all n, let Pn = \Un ' (^x\ be the progress from 
the nth ancestor of point Tq = 0. From the above analysis, the law of P„ 
under the Palm probability of converges weakly to a limit when n tends 
to infinity, and p is the mean of the limit law. A similar observations holds 
for ^1 when compared to (10): as it is the case for progress, on a long path, 
the magnitude of the hop from a point to its ancestor is "boosted" by the 
presence of its offspring. 

Note that it is also possible to derive a functional central limit theorem 
for the sequence (Tfc)^ from Theorem 17.4.4 of [19]. 

Remark 4.9. An interesting question along the line of [12] and [15] is 
whether the directed spanning forest is almost surely a tree, namely, whether 
two sequences of ancestors coalesce almost surely (in dimension 2 and 3). 
Similarly, every vertex has almost surely a finite number of successors. One 
can also expect that the directed spanning tree converges toward the Brow- 
nian web (see [13, 14, 29]). If this holds true, then the directed spanning tree 
has only one semi-infinite path, whereas in view of Theorem 2.1, the RST 
has a semi-infinite path in every direction. The edge process of successive 
ancestors in the RST converges in some sense toward the edge process of 
successive ancestors in the directed spanning stree; however, the two trees 
have a completely different topology. 



4.3. Maximal deviation. We end this section with a result on the devia- 
tion of the path from its mean. 

Let R{x) or R{X) denote the path from X = (x, 0) in the DSF with direc- 
tion —ex', R{x) may be parametrized as a piecewise linear curve {t,Y{t))t<x 
in . The maximal deviation of this curve between x' and x with x' < x is 
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defined as 
(30) 



A{x,x') := sup \Y{t)\. 

t£\x'x] 



Theorem 4.10. For all x > x' , for all e > and all integers n, 

P(A(x,x') > |X - ^ _ ^t^-ny 

This theorem is a consequence of a result on the convergence rates 

Lemma 4.11. Let {Un),n £ N, be an i.i.d. sequence of random variables 
in M.'^ and let {Xi,Yi) be the coordinates ofUi. Suppose that EXi > 0, EYi = 
and E|[/i|" < oo for all n. 

Lett>0 and Tn = J2k=i ^k- We define 



i=l 



sup 

l<k<Kt 



i=l 



Then for all e > and all integers n, 

^{Mt>t^l'^+^) = 0{t-'^). 

Proof. Let c < EXi, ¥{Kt >^) = P{E\=i^ Xi>t) = 0{t~''), from The- 
orem 3.1 of [16]. Similarly, from the same theorem, 

k 



sup 

^l<fc<[t/cj 



i=l 



> t 



1/2+e 



0(t-"). 



The lemma follows from 



sup 

l<fc<[t/cj 



i=l 



>t 



1/2+e 



□ 



Proof of Theorem 4.10. The nth ancestor of Tq = (x,0) is r„; let 
Sn = Tn — X = — J2k=i Uk and if {Xi^Yi) denotes the coordinates of C/j, let 



Kx_:c, = mil n:^Xi> 



X — X 



and Mt 



i=l 



sup 

l<k<K^ 



i=l 



We have A{x,x') = M^-x'- Since the variables {Un),n G N, are not indepen- 
dent, we cannot apply directly Lemma 4.11. Let {6k), /c G N, be the sequence 

~ Q 

given in Theorem 4.6 and Uk = X]i=0fc+i *° regenerative prop- 

erty of the sequence {9k), k £ N, {Uk),k G N* is an i.i.d. sequence; note also 
that E|C/i|" <oo. 
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Let an = sup{/c : O^+i < ra}; we have 



k=l 



i=l i=9an 



Thus, from Lemma 4.11 apphed to (Uk), G N*, we deduce the result for 
iUk),k£N. □ 



5. The radial path. For x > 0, \et X = xcx and Ro{x) or Ro{X) denote 
the path from Xq = X = [x, 0) to the origin in the RST. This path (a sample 
of which is depicted in Figure 10) may be seen as a piecewise linear curve in 
M^. We denote by H{X) the generation of Tq in the RST. As for the DSF, 
we will denote by Xq^Xi, . . . ,Xh(^x) =0 the sequence of the successive 
ancestors of Xq in the RST. Note that the points the sequence (X^ ■ e^), < 
k < H{X), is not necessarily a decreasing sequence. 




1114 



6.3 



US 



Fig. 10. An example of radial path [with its origin in (1/2,1/2)]. The initial point is of 
generation 26. 



28 



F. BACCELLI AND C. BORDENAVE 



5.1. Maximal deviation of the radial path. The aim of this section is to 
bound the tail of the distribution of the maximal deviation of the path 

(31) AW:= sup \Xk-ey\. 

0<k<H{X) 

For all t G M, we define: 

y{t) :=max(0,sup{y:(t,y) G Ro{X)}), 

z{t) := max(0, - ini{y : (t, y) G Ro{X)]), 

where the supremum over an empty set is taken to be — oo. Notice that 
y{x) =0, y{t) = for t ^ [-x,x\, suptgK2/(0 = sup{y -ey-.Y e Ro{X)] and 

A(X) = max( supy(t), sup I . 

Now let R{x) denote the path from (2;,0) in the directed spanning forest 
T^ex with direction —ex, when building this forest on the same point process 
as the RST, but with all points below the line OX removed. The path R{x) 
may be parameterized by the x coordinate of the process 

R{x) = {{t,y{t)),t<x}. 

The key observation of this subsection is the following: 

Theorem 5.1. For allt<x, y{t) <y{t). 

Proof. The proof of this result is in two steps. 

Let n~(t) = {Y £R'^:Y ■ ex <t} denote the left half-plane (x < t). We 
first prove that if the two curves share a common point T of the point 
process, namely, if {t,y(t)) = {t,y{t)) = T for some T, and if 5* G 7i~{t), 
then S ■ ey < S ■ ey, where S (resp. S) is the ancestor of T in T (resp. 'T^e^)- 

Let y = y(t) = y{t) > 0. If 5 = 5, then the property is proved. 

If 5 7^ 5" and S £ (where L_ denotes the half-plane y <0), the property 
also holds since S G (we recall that the points below OX have been 
removed by construction for T-e^). 

If 5 / 5 and 5 G Li, then 5 G ^-(i) \ B{0, |T|), and thus, S e LI D 
m{TC^{t) \ B{0, \T\)m), and the property follows in this case too since S G 

B{o,\T\)nn-{t). 

To complete the proof, we now show that Ro{x) n R{x) C N U {O} U 
{{x,0)}. Suppose that the two curves intersect in the segment between the 
vertices Ti = (xi,yi) and T{ = {xi,y'i) in T^e^ and the segment between the 
vertices Ti = {xi,yi) and T2 = (x2,y2) in We may assume xi > X2 and 
xi > X2- Moreover, if L denotes the closed half-plane below the line {Ti,T2), 
then fi G L" and 12 G L. 
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Suppose that xi < xi. Note that the ancestors of Ti (i.e., T2) and Ti 
(i.e., T2) are both in 7i~{xi). From the construction of T-e^, for any point 
Y of the point process lying in Tl{xi), T2 G B{Ti, \Ti —Y\), where ~ denotes 
the closure of a set. Using the fact that xi < xi, we deduce that B{Ti, \Ti — 
T2\)r\n-{xi)r\Lc B(Ti,\Ti-T2\)r\n-{xi)r\L. Now, by construction of 
the RST, B{Ti,\Ti-T2\)nn-ixi) nLnN = {T2}. It follows that = Tg. 

If xi>xi, then simple geometric arguments show that, necessarily, T2 = 
Ti. □ 

Corollary 5.2. For all e > and all integers n, 
F(A(X) > = 0(|X|-"). 

Proof. Let a; > 0, X = (x,0) and R{x) be the path from X in the 
directed spanning forest with direction — e^, obtained when using the same 
point process as in the RST but when removing all points above the line 
OX. The path R{x) may be parameterized as R{x) = {(t,—z(t)),t < x}. 
By symmetry, the processes z and y have the same distribution and from 
Theorem 5.1, z{t) < z{t). Therefore, it is sufficient to prove the following: 




A similar statement for the path R{x) in the directed chain built on the 
whole Poisson point process was proved in Theorem 4.10. The difference 
between R{x) and R{x) appears on the vicinity of the line OX. Thus, the 
supremum should not be much affected. 

To formalize this idea, consider the event, for 3; > 0, 7 > 0, 

A^,^ = {3Y £ with |y| < 2x and \Y - A{Y)\ > x^}. 

Following the proof of Lemma 5.2 of [17], we have 

(32) P(A^,,)<Ciexp(-C7ox2T). 

Fix < e" < e' < e and 7 < ^. On A^^^^, if suptg[„a, ,j] y{t) > then 

for X large enough, there exists at least one point of R{x) n with its y 
coordinate in the interval [x^/^"*"^ ,2x^^'^^^ ]. Let Xq = {to,yo) be the right- 
most such point. Note that |Ao| > x^/^"^^ > ^/x. The path R{x) coincides 
with the path R{Xo) on the interval [— x,Xo], provided the maximal de- 
viation of R{Xo) is less than ; indeed, the infimum of R{Xo) will 
then be lower bounded by yo — x^^'^^^ > 0. The event {R{Xq) > x^/^"*"^ } 
is a subset of = UYeNnB(o,cc)\Bio,Vl^)mYl -\Y\){R{Y)) > x^/^+^"}. 
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From Theorem 4.10 there exits a constant C (depending on e") such that 
P(A(|y|, -|y|)(i?(y)) > |y|V2+e") < C7|y|-". Now, letting Z > e, we get 

me",.) <A U {A(|y|,-|y|)(i?) > \Y\'i^+'"}\ 

\YeNnB{0,x)\B{0,^/lE) ) 
< F{N{B{0,x)) > /vrx^) + lTrx^C\x\-'^/^ . 

Thus, ¥{A^ri^ri.) = 0(|x|~") for all integer n. Using this last remark, the con- 
clusion follows from equation (32) and the inequality 

f( sup y{t)>x^^^~^") <n^'r,x)+nA",x)■ 
\te[-x,x] / □ 

For X £ N, let Ront{X) be the set of offspring of X in the RST, namely, 
the set X' £ N such that X G Ro{X'). We now state a definition introduced 
in [17]. 

Definition 5.3. Let / be a positive function on M4.. A tree is said to 
be /-straight at the origin, if for all but finitely many vertices 

Ront{X)ciC{XJ{X)), 

where for all X G and e G M+, C{X,e) = {Y £ m?:9{X,Y) < e} and 
6(X,Y) is the absolute value of the angle (in [0,7r]) between X and Y. 

Theorem 5.4. The number of points of N such that A{X) > 
is a.s. finite. The RST is f -straight at O for f{X) = |X|-V2+£ for alle>Q. 

Proof. We first prove that the number K of points T„ of N such that 
A(r„) > |r„|i/2+^ is finite. 

From Corollary 5.2, for all X G M^, P(A(X) > |X|V2+e) < 1 C\X\~'^ , 
where C depends on e. From Campbell's formula, 

Eii: = E ^ i(A(T„) > |r„|^/2+e) 

POO 

= 27r / P(A((2;, 0)) > x^/^+e)^ 

JO 

;>oo 

00. 

Jo 

The rest of the proof uses the same argument as Theorem 2.6, Lemma 
2.7 of [17] (with 3/4 replaced with 1/2). □ 

/-straight trees have a simple topology described by Proposition 2.8 of 
[17] and restated in Theorem 2.1. 
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5.2. Law of large numbers on the radial path. The aim of this section is 
to prove Theorem 2.5. We use the notation introduced at the beginning of 
Section 5. Let X = {x, 0) and Xq = X, Xi, . . . , X^(^x) = O be the successive 
ancestors of X in the RST. 

The edge between Xi and is Ui = Xi-i — Xi. Letting C/q = —Xq, we 
have 

n n 

Xn = -Y.Ui and Xn-Xk = - ^^ 

1=0 l=k+l 

We define Dq := B{Xo, \Xo\) and for n > 1, 

{n-l ^ c 

y B{Xn-Xk,\Uk+l\)\ . 
k=0 J 

The density of Un+i given {Uq, . . . ,Un) is 

where c is a normahzing constant. 

Letting g be a measurable i-^ M function, we suppose that, for ah 
Y G M^, < max(C, some positive constants C and a. We define 

the statistical average of g on Ro{X) as 

H(X) 

SHiX)= X! 9{[Xn~l - Xn]x„-i), 
n=l 

where [U]v is the vector U rotated by an angle —9 and V = rcosO ■ + 
rsin^ • Cy {[U]v is the vector U expressed in the local coordinates of V). 
We will first prove that 

where vr is defined by equation (26). 

Results on the matter will follow from the corresponding results on the 
directed path. The proof of the theorem is in three steps. 

5.2.1. Step I: pseudo-regenerative times. In this part we associate to the 
RST a sequence analog to the sequence given by Theorem 4.6 for the DSF. 

Let J^n be the cr-field generated by Xq, Ui, . . . ,Un. Clearly, Xn S J^n- We 
build a coupling on the process {Un} in the same vein as in [3] for Markov 
chains. As in the DSF, we define 

r„,+i=inf|m>r„:S(0,|X^|)n| (J B{X^ - Xk,\Uk+i\) = 0^^, 
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with To = 1 and = {0,Xr„_^ - X^-^.^+i, . . • ,^r„-i - X-r„}. For n large 
enough, = {0}. Let denote the vector in the following basis: the 
origin is in ^t„_i and the basis vectors are {exj. ,ey ), where Xt-„_-^ is 

the point of coordinates {—Xi^_j^,xil^_-^). 

We recall some notation introduced in the proof of Theorem 4.4. For 
$ = {0,yi,...,y^($)} in Mx and r(^>) > 1, we define (i?($),e($)) as the 
coordinates of i^T-($)- The set S G B{Mx) is defined by 

(35) cS = {0}U{$GA/'A'\{0}:|e($)| < 7r/6, i?($) < 1}. 

Let ^ denote the sub cj-algebra of B{J\fx) generated by the sets {t(<I>) = 
1} n {Xi($) eB},B € C, where C = {{r,e) G . ^ ^ [-7r/3, vr/S], r < 1}. 
We have the following lemma: 

Lemma 5.5. There exist positive constants A and g, and a measure v 
with support on Q, such that, for any xq large enough and for all n: 

(i) P($„+iGcS|.F,J>^A>0; 

(ii) if \X^^ \ > xo and G S, then P($„+i G ^| JVJ > g^^iA). 

Proof. Define Ca{R) := {X = (r, 6*) G : 6* G (a,0],r < i?} for a < 0, 
and Ca{R) ■■= {X = {r,e) gM.'^ -.O £ [0,a),r < R} for a > 0. We may prove 
as in Lemma 4.2 that the RST is such that either c^/g(|X„|) C -D„, or 

c_Tr/Q{\Xn\) C Dn, where Dn is the set defined by equation (33) re-expressed 
in the above referential. Hence, 

P(^„+i G c,/6(l) U c_^/Q{l)\J'n) > i(l - e-^/'') = A > 0. 

The proof of property (i) is then similar to that of (23) in Lemma 4.5 for 
the RST case. ^ 

For xq large enough, if \Xt-^ \ > xq and G S, the cone C is included in 
Dm- As in the proof of Theorem 4.4, we can check that u = PqIq satisfies 
property (ii) for a suitable choice of g, where Pq was defined in Theorem 4.4. 
□ 

Fix Xq and define 

(36) n(X) := sup{n : > Xo}. 

Following the proof of the regeneration lemma [3] , we can show the existence 
of a sequence of stopping times Ni^ (on an enlarged probability space) such 
that No = and such that, for ah A G B{Mx), 

P($„+i G A,Nk = n,Tn< n{X)) = iy{A)¥{Nk = n,Tn< n{X)). 
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From property (i) and (ii), {N^} is a renewal process and it can be built 
such that Nk+i — is a geometric random variable of parameter Xg > 0. 
Let 

K'{X) ■.= snp{k:\Xr^J>xo}. 

For 1 < A; < K'{X) + 1, ■ ■ ■ ,^N^.+i) depends on the past only through 

l^Tjsi^, ^ I - Thus, we have found some homogeneous regenerative structure on 
the nonhomogeneous process {Un)- 
Let 

(37) Ok = TN^_^ 
and 

(38) K{X)=TN,,^,y 

Note that Lemma 4.3 is still valid for the RST. Indeed, this lemma relies 
on the bounds Lj and Pi which also hold for the RST. Finally, we have the 
following decomposition: 



Lemma 5.6. For all xq large enough, there exists an increasing sequence 
of finite stopping times {6k},k (on an enlarged probability space) such 
that 9q = 0, for all k > 0, F{6k+i — 9k> n\Tg^) < Ciexp(— Con), for some 
positive constants and for <k < K{X), [Uq^j^i, . . . , Uq^^^) depends on Tq^ 
only through iXg^^l. 

Lemma 5.7. Almost surely, H{X) and K{X) tend toward infinity with 
\X\. 



Proof. For /? > 0, x > 0, let Ap^^ = {3Y e N, with |y| < 2x and \Y - 
A{Y)\ > x^}. The proof of Lemma 5.2 of [17] implies that P(A/3,x) = 0{x-'^) 
for all q. 

On A"^ , H{X) > L|^|^"^J • Let M be the number of points in N such 

that H{X) < From Campbeh's formula, EM < 2tt F{A(^^^)x dx < 

oo. Hence, M is finite almost surely. 

To prove that K{X) is finite, notice that 

FiOk+i -9k> n\Te^) < Ciexp(-Con), 

which implies that a.s. Iimsup^0fc/A; < Ci < oo. Hence, a.s. for k large 
enough, I^LCifcjl ^ implies that K{X) > k. Pick k = [|-^|^~^J and the 
proof is similar to the proof for H{X). □ 
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5.2.2. Step II: identification of the limit. For < fc < K(X), let Vn 
[Un]x,,.^ and 



S,iX)= 9{Vn). 



The next step is to identify lim|x|_^oo E,Sk{X) . Notice that in the DSF with 
direction — e^, the same construction can be done on the path starting from 

aDSF 

O; we will denote similarly C/^^^, 9^^^ and 5f = E ^.dsf , , 

These sequences are i.i.d. sequence and, since X = xex, Theorem 4.7 gives 

ES^3^ = 7:{g)e, where 6 = E(efff - Oj^^^). 

From Lemma 5.6, using the Cauchy-Schwarz inequality, we get 

CO / k \ 

ESiiX) = E E 1(^1 = k)Y^ 9{Vn{X)) 

k=0 \ n=0 / 



(39) 

< 5] Ci exp{-Con)VEg\Vn{X)). 



oo 



n=0 



Since g{Vn{X)) is stochastically dominated by max(C, Z"), where ¥{Z > 
t) < e"'^* it follows from Theorem 2.4 that 

lim Eg{VniX))=Eg{U^'^). 

\X\^oo 

Since F{\Un+i{X)\ > t\J^n) < e-^*'/i2, Eg'^{Vn{X)) is uniformly bounded. 
So, from the dominated convergence theorem and equation (39), we deduce 
that lim|jj(^|^_|_oj3 ES'i(X) = Ti{g)9. 

Similarly, the law Sk{X) depends on X and k only through |^6»(.|) which 
tends toward infinity as \X\ tends toward infinity; it follows a.s. 

(40) hm E{Sk{X)\:Fe,)=^{g)e. 

\X\^+oo 

5.2.3. Step III: convergence. We may write 

There is almost surely a finite number of edges in i?(0,xo); thus, from 
Lemma 5.7, the second term tends a.s. to 0. From equation (40), we also 
obtain that a.s. lim|x|^oo 1 A'(^) Efif ^ E(Sfc(X)|J^0j = 7r(c/)^. 
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Moreover, using the same arguments as in equation (39), we get K{\Sk{X)\'^) < 
M < oo. We can then apply Theorem VII. 9. 3 of [11] which gives 

K(X) K(X) 

For g = 1, we deduce that a.s. H{X)/K[X) tends to 6. This ends the 
proof of equation (34). 

To complete the proof of Theorem 2.5, it remains to show that if g is 

Hi 

continuous, lim|x|-»+oo '^/H{X) X]fc=i di^k-i — ^k) = '^{g)- This is a con- 
sequence of Theorem 5.4 and equation (34). Indeed, a.s. for all k, as \X\ tends 
toward infinity, — Xk]xk-i tends toward Xk-i —X^. If g is continuous, 

we deduce that l/H{X) Efif ff(^fe-i - ^fc) and l/H{X) Efif ^ g{[Xk-i - 
Xk]xi,^i) have the same limit. 

Corollary 5.8. The following a.s. limits hold: 

|X|^+oo \X\ p 
H{X)-l 

where p and are defined in Corollary 4.8. 

5.3. Shape theorem. We define Gk := \T{k)\; Gk is the size of the ball of 
center O, and radius, k for the graph-distance on the RST. 

The main aim of this section is to prove Theorem 2.2, and in particular 
the fact that Gk/k"^ a.s. tends to a constant when k tends to oo. In the 
literature this constant is known as the volume growth. 

The intuition behind Theorem 2.2 is as follows: from the results of Sec- 
tion 5.2, a point k hops away from the origin is asymptotically at Euclidean 
distance ~ kp from the origin. The disc of radius d^ contains ird^ vertices 
asymptotically. 

In order to prove Theorem 2.2, we need an estimate of the tail of H{X) 
around its mean. In Figure 11 is shown a realization of T(IOO). 

The proof of the next theorem is the heart of the proof of Theorem 2.2. 



Theorem 5.9. For all q <p, there exists positive constants Gq and Gy. 

IP 

Similarly, for q> p, 



P( if(X) > H ) < Ciexp(-C7o|X|). 



F[H{X)<^-^] <Giexp{-Go\X\). 
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The proof is again in three steps. 

Let Xq = X, Xi, . . . , Xjj^^x) = O, the successive ancestors of X in the radial 
spanning tree. For 1 < n < H{X), we define the progress 

Pn{X) = |X„_i| — \Xn\ 

and for n > H{X), Pn{X) = 0. By definition, for all integers m, 

F{HiX) > m) =^[Y.Pn{x) < \X\^ . 

In order to derive the tail bounds of the theorem, we will analyze the 
asymptotic behavior of Pn{X). 

5.3.1. Step I: pseudo regenerative times. We will use the notation of Sec- 
tion 5.2 and the sequence of stopping times = ^ij • • • given by Lemma 
5.6, which are such that, for < /c < K{X), {P0^+i{X), . . . , Pe^^^{X)) de- 
pends of ^01^ only through l^efel- 

Let 

P,{X)= PniX). 
n=6»fc+l 

Prom equation (40), 

(41) lim ^Pk{X)=pe = p/5, 

\X\^oo 
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where 9 = 6^^ = \im.\x\^oo^(^i- 

Lemma 5.10. There exists positive constants Co, Ci such that, for all k, 
nPk{X)>t\Te,) < Ciexp(-Cot). 

Proof. We will prove this result for k = 1; the proof extends to all k 
with minor changes. From Lemma 5.6, the sequence [Ok) is such that, for 
some constants Co, Ci, P(6lfc+i - 6^^ > K\re^) < Ci exp(-Co/s:) for ah K>Q. 
Let s > 0. From the Cauchy-Schwarz inequality, 

K=l 

oo J, 

(42) < ^ v/p(0i=i^)Ee2^^n=i^"W 



ii:=i 

oo 



< Ciexp(-CoK)v/Ee2^^^:=i^-W. 
We have already noticed that 



E(P„+i(X)|.F„)</x. 
Calculations based on Chernoff 's inequality lead to 

Ee^^Sti < e^^^(l + Cis)^e^^?^' < Cie^i^^, 

for all s < 1. For s sufficiently small, sC( < Cq, where Co is the constant in 

(42). For such a choice of s, we get Ee"^^'-'''"^ < Ci from equation (42). The 
lemma then follows from Chernoff 's inequality. □ 

5.3.2. Step II: case q < p. Let q < p' < p and e > be such that q < 
g^p' < p, with 6 = 9~^. From equation (41), there exists xq such that for 

\X\>xo,E{Pk{X)\\Xe,\>xo)>^^: 
F{H{X) > m) 



J2Pn{X)<\X\ 



(43) 



' m—l 



^ P E Pn{X) < \X\ - xoj + nX{B{0, Xo) > I)) 
<w(Y,Pn{X)<\X\-xo] +exp(-lln- 



. n=l 



evrx? 



38 



F. BACCELLI AND C. BORDENAVE 



where we have used the classical inequality 



F{N{B{0,xo)) > I) <exp[ -lln 



eirxf 



(Lemma 11.1.1 of [28]). 

Let n{x) and K{x) be defined as in (36) and (38), respectively. Let 
Zfc be the conditional random variable Pk{X) given K{X) > k. We have 
E(Zfc|jFg^) > Note that, for all integers j. 



1 

,n=l 



J2Pn{X)<\X\ 



Xq 



(44) 



< 



Indeed, {EkS'^'^ Pk{X) < \X\-xo} is a subset of Ui<fc<(5+e)i{^(^) ^ 
Now, let X, j such that p' j — x > aj for some a > 0: 



= P(^^P„(X) < \X\-xo,n{X)>jj 

/l{S+e)j\ \ 

E Pk{x)<\x\-xoj+n0iis+e)n<j) 

/l{S+e)j\ \ /l{S+e)j\ X 

y E ^A:<i^iJ+p(^ E 0k-0k-i<jj- 



L('5+£)iJ 



< 



IXI 



{S + e)j {5 + 



P 



I 1 LC-s+ebJ 



L(<5+^)iJ 



> 



jV - 1^1 



THE RADIAL SPANNING TREE 39 

The process {Z^ — E,{Zk\J^Nk))k satisfies the hypothesis of Lemma A.l (in 
the Appendix) with to = . We deduce that 

(45) 



Similarly, 9k — 9k-i satisfies the same light tail hypothesis as and for 
xq large enough, E(0fc - \ > xq) > to > l/{6 + e). We get 

(46) P ^k-Ok^i<j]<Cie^'^. 

^ k=i 



Now let m = [|J , I = [^^J and j = m — I. For \X\ large enough, the in- 
equality p'j — \X\ > aj holds for some positive a. Putting together equations 
(43), (44), (45) and (46), we obtain that, for \X\ large enough, 

inequality for all X. 



Since ¥{H(X) > -Li) < 1, by increasing arbitrarily Ci, we get the above 



5.3.3. Step III: case q>p. The bound F{H{X) < -1^) for g > p is ob- 



1^1 

V-'-' v^*-; ^ 

tained with the same type of decomposition: 
F{H{X) < m) 

= p(^f2^p^{x) = \x\^ 

<p(^f2^Pn{X)>\X\-Xo^ 

/[{S-e)mi \ 

E Pk{X) > \X\-xoj +nOl{S-e)m\ >m) 

/L(<5-e)mJ \ /L(5-e)mJ \ 

E Zk>\X\-xoj+F\^ J2 0k-9k-i>mj 

< (J^g~Co{\X\-xo-mp') _|_ 

for p' > p and all X, m such that |-^| — mp' > am, for some a > 0. Taking 
m = [|X|/qJ and q > p' > p,we obtain the upper bound for F{H{X) < \X\/q) 
and we are done with the proof of Theorem 5.9. 
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Note that the statement of the theorem is not tight since we do not bound 
the tail of \H(X) — \X\/p\. It seems possible to derive some concentration 
inequalities. 

Remark 5.11. Letting S{X) = -jj^J^niV fi^n), the same type of 
theorem holds for any function / which is growing to infinity slower than 
exponentially, with 1/p replaced by '7r(/)/p. 

5.3.4. Proof of Theorem 2.2. We begin with the proof of equation (2). 
Notice that 

Gk=Yl HHiT)<k). 

Letting e G (0, 1), we may write 
\Gk-N{B{0,pk))\ 

-hi^ ^ 5(0,pA;) n H{x„) < k) 

+ t{Xn G B{0,pk) n H{Xn) > k)] 

n / 

N{B{0, (1 + e)pk) \ B{0, (1 - e)pk)) 



+ 



+ ^ E e B{0, (1 - e)pk) n H{Xn) > k) 

< Ik + Jk + Lk- 
From Campbell's formula and using Theorem 5.9, 

27r 

Wk) = jo F{H{xe,,)<k)xdx 
k^ Jn+£)pk 



2tt 
27r 
27r 



k^ J{l+e)pk 



From Markov's inequality, P(4 > 0) = P(4 > l/A;^) < A;2e(4) and from the 
Borel-Cantelli lemma, we obtain that almost surely Ik is equal to for k 
large enough. 
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Similarly, let p' be such that q = p{l — e) < p' < p, chosen as in the proof 
of Theorem 5.9. Using equations (43), (44) and (45) we get similarly 

-^[j a;Cie-^«(^'f-^)xdx + Cie-^"'=y xdxj 

< 'ic,e~-o^ 

Using Markov's inequality and the Borel-Cantelli lemma, we deduce that, 
almost surely, Lk is for k large enough. 

The ergodic properties of the PPP imply that 

= k~^N{B{0, (1 + e)pk) \ B{0, (1 - e)pk)) 

converges almost surely and in mean toward Airpe. Notice that N{B{0, (1 + 
e)pk) \ B{0, (1 — e)pk)) is not an increasing sequence of convex sets. To 
prove this convergence, we need to use the independent properties of the 
PPP. 

We have proved that, for all e > 0, almost surely, 



r \Gk-N{B{0,pk))\ ^ ^ 

limsup < iirps. 

k 



Hence, almost surely, 

G,^ iV(mp))^ 

The proof for the convergence is a consequence of the dominated con- 
vergence theorem. 

Equation (1) holds since we have seen that a.s. for k large enough Ik and 
Lfc are both equal to 0, where 1^ is the cardinal of T{k) f] B{0,p + ey and 
Lfc is the cardinal of T{kY H B{0,p — e). 



6. Spatial averages of edge lengths. Consider the total edge length of 
the RST for points included in the ball B{0,x): 

J2 HX€B{0,x))\X-A{X)\. 

It is well known that for the minimal spanning tree, the subadditive ergodic 
theorem allows one to prove that ^ tends almost surely toward a constant 
(see [26]). We prove that the same holds for the RST (with, of course, a 
larger constant). 
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From Campbell's formula, 

Jo 

With the change of variable u = j^, this leads to 

r fi 

= 27r / uE{L{xu)) du. 
Jo 

The dominated convergence theorem together with equation (10) gives 

/"^ 1 , r- 

lim = 27r / u—;=du = -K y2. 
x^oo Jo V2 

We will prove a stronger result: ^ converges almost surely and in 
toward ir/\/^. 

To prove this, we consider here a slightly different problem. We sample n 
points uniformly and independently on the unit disk. This defines a finite 
point set F„ = {0,Xi, . . . We can then construct the RST associated 

to this point set. The total edge length of this RST is 

N 

C{FN):=Y.\^k-AiXk)\. 

k=l 

First notice that £ is homogeneous of order 1: for all sets Fn as above 
and all positive real numbers r, we have C{rFn) = rC{Fn), where rF„ = 
{0,rXi, . . . ,rXn}- Using this and the fact that the ratio N{B{0,x))/x'^ 
tends a.s. toward vr, it is easy to check that ^ converges a.s. toward 7r/\/2 

CiF ) j — 

if and only if tends to w^. Since we have already computed the mean 

of it is sufficient to prove that converges a.s. toward a constant. 

For proving the last property, we use the smoothness of C{Fn) (Theo- 
rem 6.1) and the Rhee and Talagrand concentration inequalities. 

Theorem 6.1. There exists a positive constant Ci such that, for all 
finite subsets F and G as defined above, 

\C{FuG)-£iF)\<Ci^\, 

in particular, C{F) < Ci^/\F\ . 

Proof. C clearly satisfies the subadditive property: for all finite subsets 
F and G, 



(47) 



C{FUG)<C{G) + C{F). 
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From Lemma 3.4.1 of [27], we deduce that there exists a constant Ci such 
that C{F) < Ciy^\F\. Subadditivity then imphes 

£{F) > C{F\JG)-C{G) 

> C{F \J G) - Ci^\. 

It remains to prove that C{F) < C{F U G) + Ciy^\, for all finite sets F 
and G as above. 

Let y € G and suppose that the points Xi, . . . , X„ of F all have Y as an 
ancestor in the RST built over the set {OjUFuG. In particular, \Xi\ > \Y\. 
Suppose \Xk\ > then 

\Xk-Y\''<\Xk-Xj\'' 

= \Xk - Y\'^ + \Xj - - 2\Xk - Y\\Xj - Y\ cosXkYXj. 

Thus, if XkyXj < f , \Xk -Y\< \Xj - Y\. The inequality \Xk - + 
\Xj - yp - 2\Xk - Y\ \Xj - Y\ cosXkYXj < \Xj - Y\^ holds for \Xk -Y\e 
[0, 2\Xj - Y\ cosXkYXj] D [0, \Xj - Y\]. It follows that 

\Xk\>\Xj\ and X^j<^ 

(48) imply 

\Xk -Y\<\Xj- y I and \Xk - Xj\ <\Xj -Y\. 

Let Q{Y,X) denote the oriented angle between OY and OX. Due to 
the origin, if X is connected to y / O in the RST, then |e(y,X)| < 7r/2. 




Fig. 12. On the left, the dashed area is Q{Xi), on the right, the set of connected points 
to Y. 
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We index the n points of F connected to Y by their increasing oriented 
angle Q{Y,X) such that |0(y, Xo)| is minimal. Xi, . . . ,Xd are the points 
counted clockwise from Xq: e{Y, Xq) < Q{Y, Xi) < ■ ■ ■ < e{Y, Xi) < n/2 
and X^i, . . . ,X_(i' are the points counted counter-clockwise from Xq with 
n = d + d' + 1 (see right picture in Figure 12). We need a tie-breaking rule: 
if Q{Y,X) = Q(Y,X'), the point with the higher norm has an index closer 
to O. 

Assume now that, for a given j > 0, we have both > \Xj\ and 

X]YXj+i < |. Then from equation (48), \Xj+i -Y\< \Xj - Y\, Xj+i be- 
longs to the set G{Xj) of points closer to Y than Xj and with a norm larger 
than Xj. However, by elementary considerations, X £ G{Xj) implies that 
|6(y, X)| < |6(y, Xj)| (see left picture in Figure 12) and this contradicts 
Q(Y,Xj^i) > @{Y,Xj) (the strict inequality comes from the tie-breaking 

rule). Similarly, for j < 0, if Xj^^^Xj < f , then \Xj^i\<\Xj\. 

There are at most 6 points such that XiYXi^i > ^ (i > 0) or XiYXi-i > 
f (i < 0). Therefore, from equation (48), there are at most 6 points such 
that (i > and < \Xi\ and - Xi\ < \Xi - Y\) or (z < and 

> \Xi\ and \Xi^i — Xi\ < \Xi — Y\) does not hold. 

Let ApiX) denote the ancestor of X in the RST built on the set {O} U F. 
We have 

C{F)=Y: \X-Af{X)\ 

= J2 HAf{X)=Afug{X))\X-Afug{X)\ 

X&F 

+ y: e \x-ahx)\ 

YGG xeFnA-laiY) 

<£(FUG)+E E \X -Af{X)\-\X -Y\. 
yeCxeFnAplai^) 

If F n ApIjq{Y) = {X_ii' , . . . , Xq, . . . , Xd}, we have seen above that for at 
most 6 points \Xi- AFuG{Xi)\ = \Xi-Y\ > \Xi-Xi_i\ > -^^(^i-i)! 
(i > 0) or \Xi-Y\ > \Xi+i-Xi\ > \Xi+i - AF{Xi+i)\ {i < 0) does not hold. 

Henceforth, if J denotes the set of points such that the preceding inequal- 
ity does not hold and H{Y) = {Xd,X^d'} U J, we have \H{Y)\ < 8 and 

E \X-AFiX)\ 

xeFnA-la{Y)\H(Y) 

d d' 
i=l 1=2 
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< \X-Afug{X)\. 

H = Uygg -f^(^)> we have \H\ < 8|G|. Using subadditivity, we deduce 

£{F)<£{FuG) + £{H)<£{FuG)+C^f8\G\. □ 

Theorem 6.1 ensures that we can apply the Rhee and Talagrand con- 
centration inequahties to the functional C (Theorem 1 of [25] and Theorem 

11.3.2 of [28]): ^[J^^ converges a.s. toward its mean. Finally, we have proved 

that a.s. and in L^, 

(49) lim ^ = tt/V2. 

More generally, we may also consider a power-weighted edge and for a > 0, 
define 

£(")(F) := 

The proof of Theorem 6.1 is unchanged if we replace by £ (however, 
the constant Ci does depend on a). Define 

(50) Aa = a/ r^-^e-""' /^dr = E\0-AeAOT. 

Jo 

From equation (10), using the Campbell formula, we get that 

lim E£(")(x)/x2 =7rA„. 

We finally deduce the following: 

Theorem 6.2. For all a>0 a.s. and in , 

lim — TT = nXa. 

X— >oo X 

We can rewrite this result as l/N{B{0,x))Y^xeN\X — A{X)\"' tends a.s. 
toward A^. That is, the spatial average of the power- weighted lengths of the 
edges tends toward the distribution of the length of (0,^ei(0)) in the DSF. 

It is crucial to notice that Xa < la- the spatial average and the average 
along a long path do not coincide. 

Remark 6.3. Theorem 2.1 of [24] gives the weak convergence 
to ttAq. We have done our analysis on the length of an edge, the same type of 
result could be obtained for other stabilizing functionals. In order to derive 
weak laws for stabilizing functionals, we can directly invoke Theorem 2.1 
of [24]. 
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7. Model extension and open problems. 

7.1. Greedy forests. The radial spanning tree lies in a large class of span- 
ning forests which are locally defined. We could extend the definition of the 
radial spanning tree over a point set N in as follows. 

Let I be a measurable function from to M+ and L be a measurable 
function from R'^ x R'^ to R+. Suppose that for all x,y,z,t£ N, {X,Y} ^ 
{Z,T}, L{X,Y) / L{Z,T), for X^Y,l{X)^ 1{Y), and for all compact K C 
R+,l-^{K)r\N is finite. Then we can define the following forest J- = (iV, E) : 
for 1{Y) < 1{X), (X, Y)GEii and only if y = a.Tginmz^^,i^z)<iix) L{X, Z). 

When 1{X) = \X\, L{X,Y) = \X -Y\ and O e iV, we define the radial 
spanning tree; if 1{X) = {X,ex) and L{X,Y) = \X — Y\, this is the directed 
spanning forest; if instead L{X,Y) is the length of a well-chosen cylinder, 
we obtain the Poisson forest of [12]. 

Let R{X) denotes the set of successive ancestors of vertex at X. The 
construction of R{X) can be thought of as greedy because at each vertex of 
R{X) a one-step optimization is performed. 

As an example, in the Appendix, Section A. 2, we discuss briefly the radial 
spanning tree obtained by choosing the L°°-norm. 

7.2. The radial spanning tree of a Voronoi cell. An interesting way to 
extend the RST is to consider two independent Poisson point processes. 
No = {TO}, the point process of cluster heads (we use here the terminology 
of sensor networks, which motivate this extension) of intensity Aq and A'^i = 
{T^}, the point process of nodes, with intensity Ai. The first point process 
tessellates the plane in Voronoi cells. We denote by Vn the Voronoi cell of 
point w.r.t. the points of Nq. Two forests can then be defined in relation 
with this tessellation: 

• The family of internal RSTs: the nth tree of this forest, 7^, is the RST 
built using the points of A'^i that are contained in Vn, with as a root. 

• The family of local RSTs: if node X belongs to Vn, one defines its ancestor 
as the point of (iVi U {r°}) n B{T^, \X - T^\) that is the closest to X. 
Notice that this ancestor does not necessarily belong to Vn- Nevertheless, 
this rule defines a forest too (see Lemma 7.1). One then defines the nth 
local RST tree Un as the tree which is the union of all the paths from 
nodes with ultimate ancestor T^. 

In what follows we concentrate on the second case which can be analyzed 
using the same type of tools as in Section 3. Figure 13 depicts a sample of 
such a forest. 

Lemma 7.1. Almost surely, there exists no node X of Nq such that the 
sequence of ancestors of X based on the local RST rule contains node X. 



THE RADIAL SPANNING TREE 



47 




Fig. 13. Local Voronoi radial spanning trees of 600 nodes uniformly and independently 
distributed in the unit square, w.r.t. 10 cluster heads, also uniformly and independently 
distributed in the unit square. The Voronoi cell of one of the cluster heads is depicted in 
red. 

Proof. Let {Yi}i>Q be the sequence of ancestors of X = Yq. li Yi = X 
for some i > 0, then necessarily the points of {Yi} belong to different Voronoi 
cells (if this were not the case, then the distance to the cluster head of the 
cell to which all points belong would be strictly decreasing, which forbids 
cycles). Then one can then rewrite {li,}i>o as 

{Yi}i>o = {Zo{l), Zo{no),Zi{l), Zi(ni), . . . , Zj{l), ...,Zj{k),.. .}, 

with Zi{l), . . . , Ziijii) G Wi for all < Z < j, where {VFj}j>o is a sequence 
of cells such that Wi / PV^+i for all / < j, n/ is a sequence of integers and 
Zj{]z) = X. Let Si denote the cluster head of Wi. Then the definition of the 
local RST implies that a.s. 

\X - So\ = \Zo{l) - So\ > \Zo{2) - 5o| > • • • > l^oM - So\ > \Zi{l) - So\. 
Since Zi{l) belongs to Wi, we have a.s. 

|Zi(l)-5o|>|Zi(l)-Si|. 
For the same reasons, for all I = 1, . . . ,j — 1, 

\Zi{l) - Si\ > \Zi{2) -Si\>---> \Zi{ni) - Si\ > - Si\ 
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and 

|Z/+i(l) - Si\ > |^i+i(l) - Si+i\, a.s. 

In addition, 

\Zj{l) -Sj\> \Zi{2) - 5/1 > • • • > \Zj{k) - Sj\ = \X- So\. 
Hence, a contradiction. □ 

Let Cn denote the total length of all edges from vertices in Vn- Let Eq 
denote the Palm probability w.r.t. A'^o- We have 

Cn = J2l{T^eVn)Lm, 

m 

where Lm is the length of the link that connects to its ancestor. Using 
the fact that r„\ £Vn iS No{B{T^,\T^-T^\)) = and the fact that Lm>u 
with u < |T^, - rO| iff Ni{B{Ti,u) n B{Tl \T^ - T^\)) = 0, we get from 
Campbell's formula that 

Eo 1(t4 G Vn)L^^ = 27TX, e-^°-^' e-^iM{r^^) dv^ r dr, 

with M{r,u) the lens defined in Section 3.1. Hence, 

(51) Eo(£o) = 27rAi e"^""'^' ( f e'^^^^^'^^'Ui^r dr. 

Jr=0 \Ju=0 J 

7.3. Open problems. The local geometry of the RST is rather well un- 
derstood. Unfortunately, the distribution of the degree of a vertex is still 
unknown. It would be appealing to compute this distribution at least in 
the DSF. In contrast with what happens in the minimal spanning tree, the 
degree is not upper bounded and so the moments of this distribution could 
be large. 

Properly scaled, the path R{X) of successive ancestors of X in the DSF 
converges weakly toward the Brownian motion. An interesting problem is to 
find a functional central limit theorem for Ro{X). Along this line, we may 
prove that the DSF converges weakly toward the Brownian web. Proving a 
weak limit for the RST is a challenging question. 

APPENDIX 

A.l. Distribution for the nth point of the PPP. We will denote by Tn 
the nth point of the Poisson point process when sorting points according to 
their distance to the origin. For a bounded Borel set A, let 

N{A) = Y,t{Tn£A). 
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We will denote the angle of T„ by and its norm by |T„| = Un- 

The sequence {4>n} is i.i.d., uniform on (0, 27r) and independent of the 
sequence {vn}- 

The random variable has a Gamma distribution of parameter (n,7rA) 
and the sequence {i^n} is a Markov chain with kernel of density 

K(y,t) = 27rAte"^"(*'-^'\ t>y. 

Indeed, the Markov property is immediate from the definition of Poisson 
point processes. In addition, by a direct martingale argument. 

Hence, 

„/ „,,2 , , .,.,2 Avr 



n + Avr 

so that conditioned on f„ = y, v'^+i snm. of and of an exponential 

random variable of parameter Att. Hence, v"^ has a Gamma distribution of 
parameter (n,7rA) indeed. In addition, the conditional density of Un+i given 
z/„ = y is 

Ut,y)=27:\te-^^'^^~y"\ t>y. 

Let 'yn,n be the density of a Gamma distribution of parameter (n,7r); 
letting Ln = \Tn — A{Tn)\, we have 

P(L„ > r) = P(i/„ > r; iV(5(r„, r) n 5(0, = 0) 

F{N{B{Tn, r) n B{0, Vt))= Q\ul = t)7n,,r(t) dt 

+00 / A/fl',/7 T-U '^-l 



„2 



M(Vt,r) 

1 ^j^2 ) 'yn,n[t)dt. 



The same type of computation can be done on the distribution of the edge 
(r„,^(r„)) or for Dn = ED{Tn). 

A.2. Radial spanning tree with Z°°-norm. In the Lqo case, if X = (t cos 0, 

tsinO), then for r < \X\oo = t max(| cos(0)|, | sin(0)|), 

M{X, r) = 2rH{r < g{X)) + (r^ + rg{X))l{r > g{X)) 

with 

g{X) = t(max(| cos0|, | sin^|) — min(| cos0|, | sin0|)) 
(see Figure 14). 

In addition, when denoting by L{X) the Lqo distance of point X to its 
parent vertex in the tree, we still have 

P(L(X) > r) = l(r < |X|oo)e-^^(^'^). 
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IsiiUlli^l^i) 



r 



Fig. 14. The L, 



'oo 



'Hens." 



A.3. Tail inequality for sums of nonnegative random variables. 

Lemma A.l. Let (Xfc)fc be a sequence of real random variables and 
{J^k}k OL filtration of this process. //E(X„+i|^„) =0 and^(\Xn+i\ >t\Tn) < 
Ci exp(— Cot), then for all to > 0, there exists positive constants Cq, Ci such 
that, for all t>tQ, 



Proof. This lemma relies on a classical computation on large devi- 
ations, we only give a sketch of the proof Let A„(A) = \D.W.{e^-^"-\!Fn-i) 
and A* (t) = sup_)^ At — A„(A), the Fenchel-Legendre transform of A„. With 
A*(0) =0, the condition P(|X„+i| >t|J"„) < Ciexp(-Cot) ensures that, for 
t / 0, A* (t) is positive and lower bounded uniformly in n by a positive Cq. 
For t > 0, A* is nondecreasing and convex (refer to [9]), hence, if t > to > 0, 
A;(t) > t^^i^ > tCo. As usual for upper bounds in large deviation, the rest 
of the proof follows from Chernoff 's inequality. □ 

A.4. Some tail inequalities in the GI / GI / oo queue. This section uses 
the classical notation of queueing theory. This leads to conflicts with some of 
the notation used elsewhere in the paper. This should not lead to difficulties, 
as only the results stated in Lemma A. 2 will actually be used. 

Let {(Tn,Tn},n £ Z, be an i.i.d. sequence of M+ x ]R+-valued random vari- 
ables representing the service times and inter-arrival times in a GI / GI /oo 
queue. For n fixed, the random variables (7„ and t„ are possibly dependent. 
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We set To = as the arrival time of customer 0; for n > 1, T„ = X]fc=o '^k is 
the arrival time of the nth customer. Let Y G be a nonnegative initial 
condition, independent of the {(T„,r„} sequence. We set W'q ' =Y, and for 
n > 1 , we define 

n— 1 n— 1 



W^^l = max y - V Tfc, max <Ti_i - V 



2<i<n 

fc=0 fc=i-l 



= max( y - (r„ - To), max ai-i - {Tn - Tj_i) 

V 2<i<n 
\Y] . 

The random variable Wn is the largest residual service time just before the 
arrival of the nth customer in the GI / GI /oo queue with initial condition 

y. 

The following additional assumptions are made: 

(i) There exist constants Co, Ci such that P((Ti > < Ci exp(— Cot) and 
P(y >t) <Ci exp(-Cot). 

(ii) The support of o"i is M+. 

(iii) P(ri>0)>0. 

The Loynes' sequence {M„} of this GI/GI/oo queue is defined by Mq = 
and 

Mn = max - ^k] , n > 1. 

This sequence is nondecreasing in n and it a.s. converges to 



(52) M = supL„i- V 

i<0 \ 



k=i~l / 



The random variable M is a.s. finite. Indeed, we can easily obtain a 
stronger assertion. Indeed, by assumption (i), for some s > 0, 



Eexp(sM) < l + ^Eexp[ s[cri_i - ^ r^j j 

KO \ \ k=i-\ ) ) 

< 1 + Eexp(g(Ji) Eexp(— siri) < cxo. 

From this and Chernoff 's inequality, we deduce that 
(53) P(Af > t) < Cie-^o*, 

for some constants Co, C\. 
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Now, we define 

(54) i/(y)=z/ = infjn>2:y-^rfc<o|. 

I k=0 J 

From time on, the initial workload does not count anymore, that is, for 
n > z^, wji^^ = max2<i<r!,(ci-i — r„ + Tj_i)+. Note that u has the same dis- 
tribution as 

v' = max< n<-l:Y -^Tk<0>. 



More generally. 



and 



k=0 k=l 



k=n 



,'ru-l 



H -Tfc- 

\k=-i^' k=-v'+l / 



+1 

have the same law, which implies that M^i and Wv have the same distri- 
bution. Since Myi < M, we have 

(55) 1 <st M. 

Note that this bound is uniform in Y . 

From assumption (iii), we may find c > and e > such that P(ri > 
c) > e. Let Bk = cl{Tk > c). Using the independence between r and Y and 
Hoeffding's inequality, 

P(. > n) < P X: 5, < y j < E exp (- i^^l^) . 

By assumption (i), y is such that P(y > t) < Ci exp(— Cot), for some positive 
constants Co, Ci, hence, 

(56) P(t/ > n) < P(y > nto) + expf - (^^^ " ^^o)^ ^ < Ci exp(-Con), 



2c2 



n 



for some positive constants Co,Ci, uniformly on the initial conditions Y. 
Let J-n be the c-field generated by the random variables Y and {{crk,Ti;),k = 

0, . . . , n — 1}. The sequence {Wn } is a {.F„}-Markov chain and the random 
variables 
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with i'{W) defined in (54) and with 2^1 = 7^ = I'iY), are {J'-"„}-stopping times. 
Using what precedes, one gets by induction that each is a.s. finite and 
that for all n, 

(57) P{i^n+i -vn> m\T,^) = P{un+i - ^n> m\wZ^) < Cie-^°™ Vm, 

(58) P{WZ\, > xl^.J = PiWlll > x\wll^) < PiM > x) Vx. 
Using (57) and a Chernoff type bound, one gets 

(59) P{vn > an) < Cie-^o'', 

for some positive constants a,Co,Ci. 
Consider now the {J^ri}-stopping time 

(60) e{Y) = 9 = mf{n > 1 : (a„_i - t„_i)+ = T^f 1}. 

Since Wl <st M, by assumption (ii) implies 

P{e <ui + 2)< P{M + a^^ - r^, - t^,+i < a^^+i - t^,+i) = 5>0. 

In the same vein, when using (59) and (58), one gets that 

P{e > n) < (1 - 5)L"/-J + P(z.L./aJ >n)< Cle-^o^ 

for some positive constants Co,Ci. 

We define the sequence = 0, 6*1 = ^(1^) and 

(61) en+i=mi{k>6n:ak-i-Tk-i = wj^^}, n > 1. 

Let = l(W^iri2 = ^.n+i - T.„+i))+. We have P(e„ = > S for aU 

n. Letting <e < 6, we have 



/ [n/a\ \ 

PiOn > (<5 - e)n/a) < P{i^in/a\ >n) + pi Ck<ni6- e)/a 

\ k=i J 

(62) 



<Cie 



-Con 



when using (59) and a Chernoff type bound. 
We have proved the following lemma: 

Lemma A. 2. Let 9n be the stopping time defined in equation (61). Under 
the foregoing probabilistic assumptions on Y , (tu) and (cn), there exists 
P > such that, for all k, 

P(0fc+i - 0fc > ra) < Cie^^o" and P(0„ > /?n) < Cie'^o", 

for some positive constants Co,Ci. 
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MATHEMATICAL NOTATION 

• A{X): the ancestor of point at X in the RST — see Section 2.2; 

• A'^{X): the ancestor of level k of point at X in the tree — see Section 2.2; 

• ^ei(^): the ancestor of point at X in the DSF with direction ei; 

• B{x,r): the open ball of radius r and center x; 

• C{X): the number of links of the tree that cross the circle of radius \X\ — 
see Section 3.5; 

• Cq: a small constant — see Section 2.2; 

• Ci: a large constant — see Section 2.2; 

• D{X): the degree of a point located at X — see Section 2.3.2; 

• Dn'- the domain (13) associated with the nth point of the path; 

• A(x): the maximal deviation of Rq{x) and the x axis — see equation (31); 

• A(x, x'): the maximal deviation of R[x) between x' and x defined in equa- 
tion (30); 

• {<I>n.}: Markov chain with state in the space of finite point processes — 
equation (17); 

• 7i: the right half-plane (x > 0); 

• Ti.^{t): the left half-plane {x <t); 

• H{X): the number of hops from point X to the origin — see Section 2.3.2; 

• L{X): the length of the link connecting a point located at X to its ancestor 
in the RST — see Section 2.3.2; 

• N: the underlying Poisson point process; 

• p: the mean stationary progress in a long path of the DSF — equation (27); 

• P{X): the progress of the link from a point at X toward the origin — see 
Section 2.3.2; 

• R{x): the path from (x,0) in the DSF — see Section 4.3; 

• Rq{X): the path from X to in the RST — see Section 2.2; 

• S: the small set of the Markov chain {$n} — equation (22); 

• aA'- return time of the Markov chain {^n} to set A — equation (18); 

• T: the RST; 

• Te^: the DSF with direction ei; 

• T{k): the subtree of the RST with vertices less than k generations away 
from the root; 

• {T„}: the points of N\ 

• T„: the Markovian times in the sequence of successive ancestors of a 
point — equation (15); 

• ^n'- the projection oiTi. \ Dn on the x axis — equation (16). 
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